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We propose a novel parameter, the anyonic topological entropy, designed to detect the error 
correcting phase of a topological memory. Unlike similar quantities such as the topological entropy, 
the anyonic topological entropy is defined using the states of the anyon occupations. As such, 
though the parameter deals with phases and phase transitions that are quantum in nature, it can 
be calculated solely from classical probability distributions. In many cases, these calculations will 
be tractable using efficient classical algorithms. The parameter therefore provides a new avenue for 
efficient studies of anyonic systems. 



I. INTRODUCTION 

Topological models that support anyonic quasiparticles 
have generated a great deal of interest, both in condensed 
matter physics [HI] and quantum information theory. 
This is due in the most part to proposals for fault-tolerant 
quantum computation [ST-fTT]. The topologically ordered 
states of these models cannot, in general, be detected 
by any local order parameter. Instead, non-local order 
parameters can be used, such as the topological entan- 
glement entropy fT2TfTSI . Though these are well suited 
to many problems, their calculation can be difficult in 
some cases. Also, some are defined specifically for pure 
states, and can run into problems when mixed states are 
considered pH) . As such, a wide toolkit of parameters 
are desirable, so that a wide variety of problems can be 
tackled. 

Here we propose a novel parameter that is explicitly 
defined to provide a witness for both error correct ability 
and the presence of topological order in models with any- 
onic quasiparticles. As such, we precede the definition of 
the parameter with a discussion of these issues and the 
conditions they place upon a system. 



II. TOPOLOGICAL QUANTUM MEMORIES 
AND ERROR CORRECTION 

For any topological model that can be used as a quan- 
tum memory, it is important to know where the bound- 
aries of the error correcting phase lie. To see how this 
may be done, let us consider the planar code [5J [7] as 
a concrete example. This is defined on the spin lattice 
of Fig. [TJa), where a spin- 1/2 particle is placed on each 
vertex and the plaquettes are bicoloured as p- and s- 
plaquettes. The following Hermitian operators are then 
defined around each of these, 

A s =n*f, b p =h<t*. a) 

These operators determine the anyonic occupation of 
their corresponding plaquettes. Any state within the 



— 1 eigenspace of an operator is said to have an anyon 
residing on the corresponding plaquette. These are so- 
called magnetic anyons m on the p-plaquettes and elec- 
tric anyons e on the s-plaquettes. A pair of e anyons is 
created on neighbouring plaquettes whenever a a z oper- 
ator is applied to their shared spin. Further applications 
of a z operators on spins forming a chain can be used to 
move the anyons around the lattice. Similarly a pairs of 
to anyons are created and moved by a x operations. 

Since all the A s and B p operators commute, they form 
the stabilizers of a stabilizer code. The anyonic vacuum, 
the subspace of states for which no anyons are present on 
any plaquette, is the corresponding stabilizer space. For 
the planar code, this space is two-dimensional, allowing 
a single qubit to be stored. These two states correspond 
to the anyonic occupations of the edges, which are not 
determined by the stabilizers. The left and right edges 
each hold either the vacuum or an e anyon, and the top 
and bottom hold either vacuum or an to. The X (Z) 
basis of the stored qubit may be chosen such that the 
| +) (| 0)) state corresponds to the vacuum on the top 
(left) edge and | — ) (| 1)) corresponds to an to (e) anyon. 

The effect of errors, which occur due to perturbations, 
thermal effects or other unwanted influences, is to move 
anyons, or to create and annihilate them in pairs. If 
the anyons are moved off of the edges of the code, log- 
ical errors occur on the stored information. With the 
conventions chosen above, a logical X is caused if an to 
anyon is moved off of the top edge, and a logical Z when 
an e is moved off of the left edge. When the errors are 
sufficiently weak, it is possible to perform a correction 
procedure to undo the logical errors caused. Such a pro- 
cedure involves measuring the anyonic occupation of each 
plaquette and, using that information, deducing the par- 
ity of the number of e anyons that were moved off of 
the top edge. If the number is odd then the net effect 
is that a single e has been moved off the edge, and so 
a logical X error is caused. If the number is even then 
there is no net anyon, since an even number annihilate 
to the vacuum, and so no logical error. Similarly, the 
parity of the number of to's moving off of the left edge 
can be determined to see whether a logical Z error has 
occurred. Once it is known what logical errors have been 
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The region A is the bottom of the annulus, B is the sides 
and C is the top. The topological entropy is then, 



FIG. 1. (a) The spin lattice on which mLxi planar code 
is defined, with s-plaquettes shown in blue and p-plaquettes 
in white. A spin-1/2 particle resides on each vertex. The 
linear size L is characterized the number of s-plaquettes along 
each side, with L = 4 in this case, (b) The regions used for 
the calculation of the anyonic topological entropy. Note that 
these partitions are of anyonic pseudospins, not of physical 
spins. The thickness of the region B must be on the order of 
its height, and all dimensions should be on the order of the 
system size. 



caused, they can be undone, and the stored information 
retrieved reliably. However, when the errors are suffi- 
ciently strong, such a correction procedure is no longer 
possible. The distribution of anyons becomes such that 
there is no reliable means to determine the nature of net 
anyon that has moved across any edge, and so whether 
logical errors have occurred. 

Successful error correction therefore requires there to 
exist an algorithm that, given the anyonic occupancies 
on the code, can determine the net anyon that has been 
moved over each edge. A more general statement of this 
condition is required for more complex encoding schemes, 
such as that in In general, successful error correc- 
tion requires there to exist an algorithm that, given the 
anyonic occupancies of a region bounded by multiple in- 
dependent edges, can determine the net anyon that has 
been moved over each edge. 



III. TOPOLOGICAL ENTROPY AND LOOP 
CORRELATIONS 

To determine whether the state of an anyonic system 
belongs to the topologically ordered phase, the topologi- 
cal entropy is usually used as an order parameter [T2|fT3]. 
This entropy is defined such that it detects whether loop 
correlations that are present in the state of a system. 
These are non-local correlations along closed loops that 
do not exist along open strings. They are present in the 
topologically ordered states of anyonic systems, and ab- 
sent in topologically trivial states. 

Let us consider the topological entropy of Levin and 
Wen [13j [14], and review its definition and relation to 
loop correlations. For this quantity an annulus shaped 
region of the system is considered and split into three. 



7 = Ia,bc — Ia,b- 



(2) 



Here Ia,b = Sa + Sb — SaB is the mutual information 
between the regions A and B, and Sa is the Von Neuman, 
entropy of the region A, etc. This entropy therefore quan- 
tifies the correlations that can be detected only when the 
entire annulus is considered, and is not detectable for any 
horseshoe shaped sub-region. Put another way, it quanti- 
fies the correlations that exist only around closed loops, 
and not open strings. 

It is intructive to consider the a stabilizer state of the 
planar code as a concrete example. Let us use tt to denote 
the set of all plaquettes with full or partial support on 
the spins enclosed by the annulus, and er to denote the 
corresponding set of vertices. The loop operators may 
then be defined, 



L a = l[A s , L„ = ]jB p . 



(3) 



The action of these operators on all spins enlosed by the 
annulus will cancel, leaving only non-trivial support on 
the annulus. We call these s- and p-type loop correla- 
tion functions, respectively. The L„ {L v ) operator will 
have eigenvalue +1 in any stabilizer state, and so an even 
number of spins around the loop will be in state | — ) 
( | 0) ) . The parity of the numbers of these states in A 
must then be equal to that in BC, but no such condition 
is present for A and B alone. These correlations there- 
fore contribute In 2 for each loop in the first term of 7 and 
nothing to the second. All other correlations are local, 
and contribute equally to each. This leads to 7 = 2 In 2 
for this state. 

Note that the L a and L n operators, as products of 
individual A s and B p , determine the net occupation of 
the a vertices and 7r plaquettes. Net vacuum or e (m) 
correponds to the +1 and —1 eigenspaces of L a (L„), 
respectively. Any state that has a definite net anyon oc- 
cupation for the enclosed region (like the stabilizer space) 
is therefore an eigenstate of these operators. It will there- 
fore have loop correlations that lead to a topological en- 
tropy of 7 = 2 In 2, meaning such states are topologically 
ordered. 

The situation becomes more complex when the net oc- 
cupancy of the enclosed region is not definite. For exam- 
ple, consider the application of single spin a x and a z to 
a state within the stabilizer space, with each error occur- 
ring with probability P. For P > this will yield a mixed 
state of multiple anyon configurations. If the loops are 
large, the total anyon occupancies of the region enclosed 
by the loop will be completely random. The probability 
of an odd or even number of | 1) 's around the loop will 
be equal, and so there will be no correlations in this be- 
tween A and BC. However, this does not mean that the 
loop correlations have disapeared, only that the means 
for their detection must be changed by modifying the 
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loop correlation function. To see why and how, consider 
measurement of the anyon configuration on an annulus 
surrounding the loop. If P is sufficiently small, and the 
annulus is sufficiently thick, an algorithm like that use for 
error correction can be used to determine the net anyon 
that has been moved over the inner edge of the annulus 
by the errors. Since this edge completely encloses this 
region, this therefore determines the net occupation of 
the enclosed region. Once its net occupation becomes 
definite, the resulting state is an eigenstate of the loop 
correlation functions and hence topologically ordered. If 
P is too high, however, the net anyon occupancy of the 
enclosed region can no longer be determined from the 
anyons of the annulus, and hence remains random even 
after the measurement. This shows that it is no longer 
topologically ordered. 

Rather than actually performing such a measurement, 
it may be incorporated into 'thickened' loop correlation 
functions with support over the entire annulus [131 115) . 
These are defined such that a state is an eigenstate if 
the total anyon occupation of the region enclosed by the 
loop, and hence by the annulus, can be determined by 
the anyon configuration on the annulus. The existence 
of such thickened loop operators means that loop corre- 
lations exist, and hence the topological entropy will take 
a non-zero value. The state is therefore topologically or- 
dered. This is true not just for the planar code, but in 
general [15] . 

This possibility gives us a novel criterion for assessing 
whether loop correlations are present for states of mod- 
els that support anyonic quasiparticles. If, by considering 
the anyon configuration on an annulus, we can determine 
the total anyon occupancy (both e and m) enclosed the 
annulus, we know that both p- and s-type loop correla- 
tion functions can be defined for which the state is an 
eigenstate. It therefore has loop correlations, a non-zero 
topological entropy and will be topologically ordered. 

If the state is an eigenstate of only one type of loop op- 
erator, and hence only the net e or m occupancy within 
the annulus can be determined but not both, then there is 
some ambiguity. It could be that the 'classically topolog- 
ically ordered' states of [T7] are present. However, it may 
also be that the state is not only an eigenstate of loop 
operators, but also of open string operators. It will there- 
fore not have loop correlations or topological order. For 
example, we can consider the spin-polarized state with 
all spins in state | 0) . This is obviously an eigenstate of 
large p-type loop operators in the planar code, since it 
corresponds to a complete lack of m anyons. However, it 
is similarly an eigenstate of open a z strings. Clearly this 
state is topological trivial. 

The reason why this ambiguity does not exist when the 
state is an eigenstate of both types of loop operators is 
because any incomplete part of a p-type loop, which will 
form an open string, will anticommute with s-type loops, 
and vice-versa. As such, states that are an eigenstate of 
open p-type strings will not be an eigenstate of s-type 
loops, and vice-versa. A state that is an eigenstate of 



both p- and s-type loops therefore is not an eigenstate of 
either type of open string operators, and hence must truly 
have loop correlations and topological order. The same 
argument holds for all other Abelian quantum double 
models [B] , and can be expected to hold for Abelian anyon 
models in general. 

IV. DEFINING THE ANYONIC TOPOLOGICAL 
ENTROPY 

Note that the criterion for loop correlations parallels 
the condition for successful error correction in topolog- 
ical memories as discussed above. The existence of the 
algorithm required for error correction directly implies 
that thickened loop operators can indeed be defined, and 
hence that the state is an eigenstate of loop correlation 
functions. To see this, consider the annulus formed by 
the region B, depicted in Fig. [ijb). This has two edges, 
one outer and one inner. The inner edge completely en- 
compasses the region A. As such, if there exists an al- 
gorithm that can determine the net anyon moved over 
each independent edge of a bounded region, it may be 
applied to the annulus to determine the net anyonic oc- 
cupation of the region A using the anyonic occupancies 
of B. A thickened loop operator can then be defined on 
the annulus for which the state is an eigenstate. 

These properties of error correction and loop correla- 
tions allow a quantity to be defined, whose value signals 
whether the former is possible and the latter is present. 
To define this quantity, first note that any state of the 
underlying spin lattice can be completely described by 
specifying the anyonic occupancies of all plaquettes, as 
well as the internal state of the quantum memory |23) . 
Rather than considering the physical spins of the sys- 
tem we may therefore map the model to one of anyonic 
pseudospins. One of these is assigned to each plaquette, 
with the states | 0) and 1 1) used to denote that the cor- 
responding plaquette is occupied by the vacuum or an 
anyon, respectively. A final pseudospin, representing the 
state of the quantum memory, can be defined in terms of 
the occupation of the left or top edges. However this will 
not enter into the arguments below. 

Consider a partition of the anyonic pseudospins ac- 
cording to the annulus of Fig. [I] (b). Using this, ps can 
be used to denote the reduced density matrix of all any- 
onic pseudospins within the region B, obtained from the 
state of the entire system by tracing out all other pseu- 
dospins. Similarly, p^ B can be used to denote the joint 
density matrix for the anyonic pseudospins of region B 
and the net anyon occupation of region A. This is ob- 
tained from the state of the entire system by tracing out 
all pseudospins not in A or B, and then also tracing out 
all degrees of freedom within A except that of the net oc- 
cupation. The entropies corresponding to these density 
matrices, S B = ~Pb log pb and S^ B = -Pab^°&Pab-> 
can then be used to define the following quantity, 

T' = S AB -S B . (4) 
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This quantifies the information concerning the net oc- 
cupation of A that cannot be deduced from the anyonic 
occupation of B. In the error correcting phase and as 
L — > oo, knowledge of B allows the net occupancy of A 
to be deduced completely. As such, F = 0. Since this 
value signifies the ability to completely deduce the anyon 
configuration within the annulus it is a witness, and a 
definite smoking gun, for the presence of topological or- 
der. 

The uncorrectable phase is that for which logical X 
and Z errors both occur completely randomly and there 
is no way to use the anyon configuration of the bulk to 
deduce whether they have occurred or not. The former 
implies that the net anyon moved over any edge is com- 
pletely random. The net anyon occupation of the region 
A is therefore a random mixture of the four possibilities 
(vacuum, e only, m only and both e and m), and so has 
an entropy of S^ B — 2 log 2. The latter implies that 
the net occupation cannot be deduced from the anyon 
configuration of B, and so Sj g = + Sb- As such 
F = S A = 2 log 2 for this case. 

When F takes an intermediate value between these 
two extremes, we find the ambiguous case mentioned in 



Section III In the thermodynamic limit, and assuming 
local interactions, the only possible intermediate value 
will be F = log 2, corresponding to unambiguous deter- 
mination of the net occupation of only one anyon type. 
This may or may not signal that one type of loop cor- 
relation is present, depending on whether the underlying 
spin state is an eigenstate of open string operators as 
well as closed loops. In order to determine which is the 
case, additional insight into the model considered must 
be used. Otherwise, this intermediate value is neither a 
witness that topological order is present, or proof that it 
is not. 

The quantity F is atypical in that it assigns the value 
of zero to states that are ordered and a value of 2 log 2 
to those that are not. As such let us define and use the 
following modified version, 



T = 21og2-F = 21og2 + S s - Si 



AB ' 



(5) 



For this, it is the value T = 2 log 2 that is the witness of 
states in the topologically ordered and error correcting 
regime (corresponding to F = 0), and T = for those 
that are not (corresponding to F = 2 log 2). It is this 
quantity that we call the anyonic topological entropy. 

It is important to note that, though T = implies 
error correction is not possible, this only holds for the 
particular definitions of the anyons used in the calcu- 
lations. There may still be a definition of anyons for 
which r = 2 log 2, and so error correction can be per- 
formed. For example, a state stabilizer by all the A s and 
B p operators of the planar code is clearly topologically 
ordered. When the anyonic occupations are defined by 
the eigenspaces of these operators, this state corresponds 
to the anyonic vacuum and gives T = 2 log 2. However, 
if the occupations were instead defined by rotated oper- 
ators A' s and B' p , for which each a x is replaced with a a z 



and vice-versa, the distribution of anyons would appear 
completely random and uncorrectable, giving T = 0. 

Note that these arguments hold when the distribution 
of anyons is classical, but could run into complications 
for a quantum distribution (since, for example, it is pos- 
sible for Sb to be greater than Sx B ). However, since 
error correction depends only on correlations in the any- 
onic basis, calculations can be made for quantum states 
by removing the off-diagonal elements of the density ma- 
trices in this basis. The quantum states thus become 
classical, and the above may be applied normally. 

Since the definition of the anyonic topological entropy 
is based on classical distributions of anyon occupancies, 
ignoring the coherence of the underlying spins, it gives 
a significant practical advantage over other methods in 
some cases. It allows the calculations of T, or al least 
bounds on the value, to be computed using efficient clas- 
sical algorithms, such as those already used for error cor- 
rection H8H221. 



A. Relation to error correction 

The value of T is discontinuous at phase transitions 
between topologically ordered and non-topologically or- 
dered states. This discontinuity can be seen in the exist- 
ing results of error correction algorithms 19-22J. Though 
these do not consider an annulus, they equivalently must 
determine the net anyon moved off independent edges of 
a bounded region. In these the probability of successful 
correction is found to be a step function in the thermody- 
namic limit, shifting from a value of 1 below a threshold 
to 1/2 above. This then becomes a step function between 
2 log 2 and for F This is considered in more detail in 
the Section |Vl 

Since the definition of the anyonic topological entropy 
is based on error correction in topological memories, one 
might conclude that it is equivalent to a direct analysis 
of the fidelity of a quantum memory (achieved by per- 
forming error correction and determining with what rate 
this is successful in correcting errors). However, though 
they are indeed equivalent in all cases to which both may 
be applied, they have a very distinct difference which al- 
lows anyonic topological entropy to be applied to a wider 
variety of problems. The fidelity of a memory is a pro- 
cess based quantity: it only makes sense if the memory 
is prepared in an initial state to which an error model is 
then applied, such that the final state after error correc- 
tion can be compared to the initial state. The anyonic 
topological entropy, however, is purely state based in its 
definition. It can be applied to any state, whether it is 
the result of a dynamic process (such as that considered 
by the fidelity), or whether it represents an equilibrium 
for which there is no concept of an error model or initial 
and final states. The anyonic topological entropy can 
then, for example, be applied to thermal states to de- 
termine whether a memory will be stable against finite 
temperature. An example of this is considered in section 
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CALCULATING T USING ERROR 
CORRECTION 



B. Relation to topological entropy 



The anyonic topological entropy, like the topological 
entropy, is a witness to the non-local loop correlations of 
topologically ordered states. In this sense they are equiv- 
alent. As described in Sec. Ill if the net e (m) anyon 



occupation of A can be determined by the anyon configu- 
ration of B, then a modified p-type (s-type) loop operator 
can be defined for which the state is an eigenstate. This 
shows that the loop correlations exist, with each giving 
a contribution of log 2 to the topological entropy 7 (here 
we use the definition of [13] )• Though this implies that 
r = 7 (at least when the correct definition for anyons 
is used for the former), this relation must be used with 
care. The ambiguities described above mean that the 
only definite case is for T = 2 log 2. Here we know that 
both types of loop correlations are present are definitely 
present, and so 7 = 2 log 2. Otherwise T and 7 only take 
the same value when ambiguities are not present. 



Generalization to other Abelian models 



As discussed above, given knowledge of the error model 
applied to an anyonic vacuum state, and given the anyon 
configuration within the bulk, an error correction algo- 
rithm can determine the most likely net anyon moved 
off each edge. If applied instead on the annulus of Fig. 
[ljb), this means that the algorithm can determine the 
most likely net occupation of the region A when given 
the anyon configuration of B. 

Using such an error correction algorithm, we can use 
numerical simulations to calculate T. Here we demon- 
strate how this may be done for the planar code subject 
independent bit and phase errors, where the e and m 
anyons can be considered separately. Corresponding cal- 
culations for other error models and topological codes 
can be done in a corresponding manner. 

Let us use P s to denote the probability that the algo- 
rithm correctly guesses the net e occupation of A, and n s 
to denote the probability that the net e occupation of A 
is the vacuum. Also we use P' s = 1 — P s and n' s = 1 — tt s . 
These values can be obtained numerically by running the 
error correction procedure for a large number of samples 
[T^I22j . In terms of these quantities, the mutual infor- 
mation between the net e occupation of A and the guess 
of B can be expressed, 



A.B 



s(p s ir s + py s )-s(p s ). 



(7) 



The generalization of T to other Abelian quantum dou- 
ble models is straightforward. One simply replaces the 
two level pseudo-spins used in the case of the planar code 
with ones of sufficient dimensionality to record the oc- 
cupancies for the quasiparticles of the model. For ex- 
ample, for a quantum double model D(G) based on an 
Abelian group G, there are \G\ different anyon types (in- 
cluding the vacuum) that can reside on each plaquette 
[BJ. Pseudo-spins of dimension \G\ are therefore required 
to record this information. The partitions and mutual in- 
formations are then applied exactly as above. In general 
for Abelian models, the total number of anyon types is 
given by D 2 , where D is the so-called quantum dimension 
of the model. The definition of V then becomes, 



Here S(x) = — x log x — (1 — x) log(l — x) denotes a Shan- 
non entropy. Note that this value for P~ actually yields 
a lower bound, because error correction algorithms are 
not perfect. They cannot, in general, deduce as much 
information about the net occupation of A as would be 
possible via brute force methods. A similar lower bound 
1 7 can be calculated for the mutual information be- 

A,B 

tween the net m occupation of A and the guess of B. 
The lower bound for the total mutual information I 
is the sum of these two contributions. 

Since V may be expressed V = S A 
quantity may then be used to lower bound T 



A.B 



J| B , the above 



r > 2 log 2 + F Ab + P AB - S(k s ) - S(tt p ). (8) 



T = 2logD + S B -S A 



AB- 



(6) 



The value of T — 2 log D is again the witness of topo- 
logical order, since only topologically ordered states can 
realize the full confinement of anyons required for error 
correction. This value implies that the topological en- 
tropy will also be 7 = 2\ogD. T = again signifies that 
topological order is not present, and intermediate values 
are ambiguous. Generalization of the quantity to models 
with non- Abelian anyons or non-anyonic topological de- 
fects is not straightforward, but can be expected to follow 
from the same basic principles. 



This allows studies to be made of how T changes with 
error rate, or with the time spent coupled to a thermal 
bath, that correspond to the results already obtained for 
error correction [TM2"2"| . In these works, P s and P p ap- 
proach a step function as L — ¥ 00, with P sp — > 1 in 
the error correcting regime and P s ^ p — ¥ 1/2 in the uncor- 
rectable regime. The same would be true for brute force 
methods, with the only difference being that the imper- 
fect algorithms have lower thresholds and take longer to 
reach their asymptotic values j7]. From Eq.'s Q and ^ 
we see that this will lead also to a step function in T, 
with r — 2 In 2 in the error correcting regime and T = 
in the uncorrectable regime. 
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VI. APPLICATION TO THERMAL STATES 

The topological entanglement entropy of thermal 
states of the toric code have been studied previously, 
though not without somewhat complex calculations |17[ 
[25] . Here we demonstrate that corresponding studies of 
the planar code using the anyonic topological entropy are 
much simpler. 

A study of thermal states first requires a Hamiltonian. 
The Hamiltonian H of the planar code, and its corre- 
sponding form for the anyonic pseudo-spins, is, 

h = -j s J2 a *~ JpJ2 b p ( 9 ) 

s p 

hp = -jsJ2 < - J p E a l ■ ( 10 ) 

8 p 

Given this Hamiltonian, the thermal state is the tensor 
product of thermal states for each individual pseudospin, 



tr(e 



-JsP<y 



tile 



(11) 



where ft is the inverse temperature. The probability that 
an anyon exists on any plaquette can then be determined 
from these distributions to be p s = (1 +exp[2J s ft])~ 1 and 
p p = (l+exp[2J P /3])- 1 . 

Due to the fact that each plaquette's anyonic occupa- 
tion is uncorrelated to any other, it is easy to see that 
I A b = 0. As such r = 21og2 - S A = 21og2 - S{ir s ) - 
S(irp), where its (tp) is the probability of an even parity 
of anyons, and hence a net occupation of the vacuum, in 
the s-plaquettes (p-plaquettes) of region A. These can 
be easily found from the probabilities p s and p p to be 
TT S = (1 + [1 - 2p,] n «)/2 and n P = (1 + [1 - 2p p ]»*)/2, 
where n s and n p are the number of s and p plaquettes in 
the region A, respectively. 

Let us characterize the size of a code by its total num- 
ber of plaquettes, N. When calculating V for systems 
of differently sized codes, the region A should be cho- 
sen to grow with the system size. Let us choose to 
include half of all plaquettes within the region A, and 
hence n s = n p = N/4. The corresponding graphs of T 
against inverse temperature ft are shown in Fig. [2] for 
Js = Jp = 1. These graphs show remarkable agreement 
to those in previous studies for the topological entangle- 
ment entropy |171 123) . In the thermodynamic limit and 
for any finite temperature it is easy to see that the any- 
onic occupation of A will be completely random, giving 
r = 0. This reproduces the known result than the ther- 
mal state of the planar code is not topological ordered in 
this case [T71I23]. 




FIG. 2. A graph of F against temperature T = 1/(3 for 
planar codes of various sizes N. Logarithms are taken base 2. 
Topological order is found to be present for low temperature, 
but gets weaker for larger systems. 



other cases are beyond the scope of this work. How- 
ever we conclude with some discussion of how general- 
ization to other cases may be possible. The definition 
of r can be straightforwardly generalized to models for 
which point-like quasiparticles exist in a space of any di- 
mension, rather than just two dimensions as above. The 
only requirements are than the region A must be closed, 
the region B must surround A, and the composite re- 
gion AB may not cover the entire model. The latter is 
required such that no global features (such as a conser- 
vation law) can be applied to deduce the net occupation 
of A from B even when error correction is not possible. 

Notably, the parameter can be applied to cases for 
which the topological entanglement entropy cannot. The 
prime example is that of the planar code subject to peri- 
odic application of noisy stabilizer measurements, which 
gives a three dimensional (2 space and one time) configu- 
ration of defects [7112] ■ This model is the most realistic for 
quantum computation and has the highest known noise 
threshold, of around 1% 0, and so the anyonic topologi- 
cal entropy may be useful for illuminating studies of this 
important model. 

Note that in other cases for spacial dimensions other 
than two, the phases detected by the parameter will not 
necessarily be topological in nature. This is because the 
above arguments concerning loop correlations may not 
apply. An example of the application of the parameter 
to a one-dimensional system is shown in the next section, 
where it detects the boundary of a ferromagnetic phase. 



VII. POSSIBILITIES FOR GENERALIZATION 



VIII. APPLICATION TO A 
ONE-DIMENSIONAL MODEL 



The anyonic topological entropy has been defined here 
for use with Abelian anyon models, and claims about 



As an example of the calculation of T for a one- 
dimensional model, we consider the ID transverse field 
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Ising model [23], which is equivalent to Kitaev's Majo- 
rana chain [25H27] . The Ising model with no field can be 
considered to be a one-dimensional variant of the planar 
code, where the lattice used is a line rather than a square 
lattice as usual. With the field it is therefore a perturbed 
ID planar code. The Hamiltonian of the model may be 
expressed, 



H, 



N-l 



N 



(12) 



For < fi < 1/2 it is known that the chain is in a fer- 
romagnetic phase, which corresponds also to the ordered 
phase of the Majorana chain. For the case of fi = 0, 
the excitations of the Hamiltonian can be understood in 
terms of quasiparticles, which correspond to domain walls 
between regions of aligned spins. These can be used to 
calculate T. For fi > these may not be the optimal 
choice of quasipartices, and so though T > is a witness 
of order, r — is not evidence of its absence. This may 
therefore lead to us detecting the phase transition at a 
slightly lower value of /i than the true critical point. 

To calculate T, we must first define a set of pseudo- 
spins whose states correspond to the occupations of the 
domain wall quasiparticles. N — 1 such pseudo-spins are 
required, where the state of the jth pseudospin is de- 
fined as | 0) if there is no domain wall between the jth 
and j + 1th spins and 1 1) if there is. The states of these 
pseudo-spins contain all information about the relative 
orientations of the Ising spins, but not of the actual ori- 
entations themselves. To add this information, an iVth 
pseudo-spin is introduced whose state is defined to be 
the same as that for the Nth spin ( | 0) for spin up in the 
x-direction, 1 1) for spin down). The Hamiltonian Hi can 
now be rewritten in the pseudo-spin space as follows, 



N-l 



2 



N-2 

'£ i \C CC i ^ *C 



With the Hamiltonian defined, the ground state may be 
calculated numerically for different values of /i. The off 
diagonal terms may then be removed from the state, ren- 
dering it a classical distribution of domain walls, and the 
corresponding values of Y can be determined. The re- 
sults of such a study are shown in Fig. [3j It is found 



that r = log 2 at /j, — 0, showing that this state is or- 
dered as expected, and r = for fj, >• 1/2, showing that 
the order has been destroyed. Results for the system sizes 
considered intersect at around /i = 1/2, suggesting that 
this is the transition point between the ordered and dis- 
ordered phases, in agreement with known results. The 
intersection point appears to be just below the known 
critical value, which may be due to our use of fi = 
quasiparticles throughout. 

It should be noted that a similar study could be per- 
formed for the toric code perturbed by magnetic fields. 
Though the phase transitions for this case have already 
been studied [25]: the use of the anyonic topological 



o 
o 




U 








FIG. 3. A graph of Y against fi for TV = 5, N = 10 and 
N — 15. To calculate Y the chain of TV anyonic pseudo-spins 
was cut into five equal fifths. The first and fifth of these were 
used for region C, the second and fourth for B and the third 
for A. Y here is calculated using logarithms of base 2. 



would provide an independent verification of these results 
and may shed new light on the properties of the model. 
However, this would require a significant increase of the 
computation cost compared to the one-dimensional case, 
and so such a study is deferred to future work. 
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